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On the abundance of gravitational arcs produced by 
submillimeter galaxies at radio and submm wavelengths 
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Abstract. We predict the abundance of giant gravitational arcs produced by submillimeter galaxies (SMGs) lensed by fore- 
ground galaxy clusters, both at radio and submm wavelengths. The galaxy cluster population is modeled in a realistic way 
with the use of semi-analytic merger trees, while the density profiles of individual deflectors take into account ellipticity and 
substructures. The adopted typical size of the radio and submm emitting regions of SMGs is based on current radio/CO obser- 
vations and the FIR-radio correlation. The source redshift distribution has been modeled using three difi'erent functions (based 
on spectroscopic/photometric redshift measurements and a simple evolutionary model) to quantify the effect of a high redshift 
tail on the number of arcs. The source number counts are compatible with currently available observations, and were suitably 
distorted to take into account the tensing magnification bias. We present tables and plots for the numbers of radio and submm 
arcs produced by SMGs as a function of surface brightness, useful for the planning of future surveys aimed at arc statistics 
studies. They show that e.g., the detection of several hundred submm arcs on the whole sky with a signal-to-noise ratio of at 
least 5 requires a sensitivity of 1 mjy arcsec"^ at 850 jim. Approximately the same number of radio arcs should be detected 
with the same signal-to-noise ratio with a surface brightness threshold of 20 yujy arcsec"^ at 1.4 GHz. Comparisons of these 
results with previous work found in the literature are also discussed. 



1. Introduction 

The effect of gravitational lensing constitutes a unique research 
tool in many astrophysical fields, since it allows one to investi- 
gate the structures of both the lenses (e.g., galaxies and galaxy 
clusters) and the lensed background sources, as well as to probe 
the three-dimensional mass distribution of the Universe. In 
particular, one of the most spectacular phenomena associated 
with the gravitational deflection of light are the giant arcs ob- 
served in galaxy clusters, which are caused by extended back- 
ground sources lying in the regions where the lensing mag- 
nification produced by the cluster is strongest. The effective 
source magnification can easily reach ~ 30 in these cases, pro- 
viding the opportunity to detect and spatially resolve the mor- 
phologies and internal dynamics of high redshift background 
sources at a level of detail far greater than otherwise possi- 
ble. An illustrat i on of t his powerful technique was presented in 
Swinbank et alj (120071) . where a magnification factor of 16 by 
the cluster RCS 0224-002 aUowed them to study the star for- 
mation activity, mass and feedback processes of a Lyman break 



galaxy at z ~ 5, something that (without the help of lensing) 
would not be possible beyond z ~ 2 with current instruments. 

A particularly interesting application of gravitational 
lensing is the so-called arc statistics, i.e. the study of the 
abundance of large tangential arcs in galaxy clusters. Among 
other things, this abundance is sensitive to the cluster mass 
function, the cluster dynamical activity (e.g., infall of matter, 
mergers, etc.) and internal structure of host dark-matter 
halos (e.g., triaxiality and the concentration of the density 
profile), which makes arc statistics a unique tool to study 
the cluster population. Arc statistics studies at optical and 
near-infrared wavelengths have been numerous in the past 



decade on both the observational dLuppino et al. 1 19991: 



Wambsganss et al.ll2005l: iFedeU et al ll2008l) . 



Gladders et al.ll2003l: IZaritskv & Gonzalez! l2003b and theoret 



ical si des ( Bartelmann & WeissI 1994 ;_ Wambsganss et al 
19981 [Bartelmann et al.l 120031 iMeneghetti et al.l 12003 



All the currently known giant arcs come from detections 
in the optical. However, the fraction of lensed sources ob- 
served in the mm/submm wavebands is expected to be much 
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served in the mm/subm m wavebands is exp ected to be much 
larger than in the optical (IB lain Il996lll997l) : Due to the spec- 
tral shape of the thermal dust emission, the observed submm 
flux density of dusty galaxies with a given luminosity remains 
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approximately constant in the redshift range 1 < z < 8 in- 
stead of declining with increasin g distance (usually referred to 



as "negative K-c orrection", see iBlain & Longainll993 , Il996t 



Blain et alj|2002l) . This effect, together with the steep slope of 
the observed submm number counts, produces a strong magni- 
fication bias that makes submm galaxies (hereafter SMGs) an 
ideal source population for the production of lensed arcs. 

The S MGs were first detected about a dec a de ago with 
SCUBaQ (Smail et al.'|l997|; iHughes etal]ll998l: iBarger et al ' 



1998.: .Bales et al... 19991) 7 The current observational evidence 
indicates that these objects are high-redshift dust obscured 
galaxies, in which the rest frame FIR peak of emission is 
observed in the submm band. Their FIR luminosities, in the 
range 10" - IO'^Lq, are ~ 100 times higher than what is ob- 
served in local spirals. Their energy output seems to be dom- 
inated by star formation processes induced by galaxy inter- 
actions/mergers, although a good fraction (~ 30 - 50%) of 
SM Gs also seems to host an Active Galactic Nuc leus (AGN, 
e.g. lAlexander et alj2005[lMichalowski et al.l2009l) . The avail- 



I might be the progenitors 



of massive local elUpticals (iLillv et al. 


1999; 


Smail et al."2002'. 


2004 


Webb et all 20031 Genzel et al 


I2OO3 


; [Alexander et al.i 


2003 


2005t Swinbank et all 20061 20081 Michalowski et al.l 


2009 





In the ~ 30 clusters ob served with SCUBA (iSmail et al 



2002 : IChapman et al.ll2002t ICowie et al.ll2002l: iKnudsen et al 
20081), only 4 multiply imaged SMGs have been reported 



to date (IBorvs et al.l l2004t iKneib et al.1 120041; IKnudsen et al 



20081) . This extremely low detection rate is due to three ma- 
jor limitations of current 850 jum surveys: (i) very small sky 
coverage (~ 3 square degrees, including all cluster and blank 
field surveys), (ii) confusion Umited maps at ~ 2 mJy (which 
means that only the brightest members are detected) and (iii) 
insufficient resolution (~ 15") to resolve extended lensed struc- 
tures. Current efforts to increase the surveyed area at 850 /im 
include the SASsH the SCLsH and the all sky survey that will 
be carried out with the HFI bolometer on board of the Planck 
satellit^^, but their resolutions will still be insufficient to iden- 
tify lensed arcs. The only instrument that can currently provide 
sub-arcsecond resolution in submm (at 890 //m) is the SMaH 

However, the tight correlation between radio syn- 
chrotron and FIR erpission observed in star-forming galaxies 
(Ivan der Kruitlll973l:lHelou et aljl985b . provides an alternative 
way to obtain high resolution images of SMGs. Commonly 
referred to as "the FIR/radio correlation", it covers about 



five o rders of magni tude in luminosity (|Condonlll992t iGarrettI 
2002h out to z ~ 3 (iKovacs et al.ll20o4 IVlahakis et al.l l2007t 



IbaretalJl2008 : iMichalowski et al 



20091) . The standard model 



about its nature considers that both emissions are caused by 



' S ubmillimeter Common User Bolometer Array ^Holland et al.l 
1 19991) . which used to be mounted at the James Clerk Maxwell 
Telescope (JCMT) located in Hawaii 

^ SCUBA-2 All Sky Survey, a ~ 2 x lO'* square degree survey with 
a 5cr depth of 150 mJy and 15" resolution 

^ SCUBA-2 Cosmology Legacy Survey, a ~ 35 square degree sur- 
vey with a 5cr depth of 3.5 mJy and 15" resolution 

* 5cr depth of ~ 350 mJy and ~ 50" resolution 

^ The Submillimeter Array in Hawaii 



massive star formation: While young massive stars produce UV 
radiation that is re-emitted in the FIR by the surrounding dust, 
old massive stars explode as supernovae, producing electrons 
that are accelerated by the galactic magnetic field and gener- 
ate the observed radio sy nchrotron emission (iHarwit & Pacini 



19751: iHelou et al.lll985l) . Therefore, given the possible com- 
mon physical origin of both emissions, radio interferometric 
observations can be used as a high-resolution proxy for the rest- 
frame FIR emission of high-z galaxies observed in the submm. 

The advent of ALMA0 will open a new window into 
mm/submm astronomy at sub-arcsecond resolution and sub- 
mJy sensitivity, allowing the detection of resolved gravitational 
arcs produced by SMGs. Although its small instantaneous field 
of view (FOV) severely limits ALMA's survey capability, a 25 
meter submm telescope (CCAlfl) is going to be built on a high 
peak in the Atacama region to provide wide field images (~ 400 
arcmin^) with a resolution of ~ 3.5" at 350 //m. With the com- 
bined capabilities of both instruments, arc statistics studies in 
the submm might be possible. 

At the same time, radio interferometry is also experienc- 
ing major technological improvements. In particular, the VLA0 
and MERLIIvfl are cuiTently undergoing major upgrades which 
will boost their sensitivities by factors of 10 - 30 and dramat- 
ically improve their mapping capabilities. The new versions of 
these arrays (e-MERLIN and EVLA) will be fully operational 
in 2010 and 2012 respectively. In order to assess the prospects 
for the study of gravitationally lensed arcs at submm and radio 
wavelengths, in this work we report detailed theoretical predic- 
tions about the abundance of arcs produced by the SMG popu- 
lation at 850 fim and 1.4 GHz. 

The paper is organized as follows. In Section |2] we in- 
troduce all the relevant quantities that are necessary to derive 
the total number of arcs detectable on the sky. In Section [3] 
we present and discuss all the observational information about 
SMGs that is required for the subsequent strong lensing analy- 
sis: morphology, redshift distribution and number counts. A de- 
scription of our cluster population model and the way in which 
the abundance of large arcs is computed is given in Section 
m The derived arc redshift distributions and number of arcs 
are presented in Section |5] while in Section |6] we discuss how 
these relate to previous findings in the literature. A summary 
and conclusions are presented in Section |7] 

The adopted cosmology corresponds to the standard 
ACDM model with cosmological parameters inferred from the 
WMAP-5 data release in conjunction wit h type la supernovae 



and baryon acoustic oscillation datasets (iDunklev et al.ll2009 
Komatsu etal]|2009b . namely D.^^ = 0.279, Qa.o = 0.721, 
0-8 = 0.817 and Hq = /dOOkm s"' Mpc"', with /; = 0.701. 

2. Strong lensing statistics 

The choice of the best parameters to be used in order to char- 
acterize the morphological properties of long and thin gravita- 
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tional arcs is still a matter of debate. In this work, we adopted 
the quite popular choice of the length-to-width radio d, which 
has to be larger than a certain threshold c/q (usually 7.5 or 10) 
in order to consider an object as a giant arc. 

Given a set of background sources placed at redshift z&, the 
efficiency of the galaxy cluster population to produce arcs with 
length-to-width ratio d > do, is parametrized by the optical 
depth 



0.6 



4kDIJo Jo 



(M,z)crd„(M,z)dMdz, (1) 



where is the angular diameter distance to the source redshift 
and n{M, z) is the total number of clusters present in the unit 
redshift around z with mass in the unit interval around M. The 
cross section o-j^(M,z) is the area of the region on the source 
plane where a source has to lie in order to produce (at least) one 
gravitational arc with d > do, for a single cluster with mass M 
at redshift z. This depends in general on the cluster structure, 
the source properties, and the redshifts of both the cluster and 
the source. Since in realistic situations sources are distributed 
at diff'erent redshifts, we can calculate the average optical depth 
by integrating T^^izs) over the source redshift distribution piz^), 



Tdo(Zs) 



Jo 



(2) 



In this way, the total number of arcs with leng th-to-width 
ratio d > do can be calculated as (IBartelmann et al.lll998.) 



(3) 



where is the observed surface density of sources, and f,;„ = 
tj„(-i-oo) is the total average optical depth (i.e., the average op- 
tical depth with the integral extending to all possible source 
redshifts). Therefore, the number of arcs produced by a given 
population of background sources can be calculated by provid- 
ing the following observational constraints: (i) the characteris- 
tic source shape and size, which is necessary to compute the 
cluster cross sections cr^„(M, z), (ii) the source redshift distri- 
bution p{Zs), which is necessary to evaluate the optical depth 
fdo, and (iii) the cumulative source number counts A^. 



3. Characteristics of the SMG population 

3.1. Source shape and size 

In strong lensing statistics studies, it is customary to character- 
ize the size of elliptical background sources using the equiv- 
alent effective radius Rg = 'sfab, which is the radius of a cir- 
cle that has the same area of the elliptical source with semi- 
major axis a and semi-minor axis b. In addition, the orien- 
tation of sources is randomly chosen, and to account for the 
different source shapes the value of the axis ratio h/a is con- 
sidered to vary within a certain interval. The typical values 
of these parameters used in optical ray-tracing simulations are 
Re - 0.5" and b /a randomly varyi ng in the interval [0.5, 1] 
dMeneghetti et al.lboOO. 2003. 2005;,Fe"deh et al.ll2006l) . 
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Fig. 1. Effective radii and axis ratios derived for the 1 .4 GHz 
radio c ounterparts of the 12 SMGs presented in Biggs & IvisonI 
(l2008h . The red points correspond to the 890 yum counterpart 
of GN20 detected with the SMA, whose size was derived by 
fitting a Gaussian (triangle) and an elliptical disk (square) to 
the data jYounger et al .l2008h . Error bars are computed through 
standard error propagation. 

Due to the low resolution of submm single-dish observa- 
tions (e.g., FWHM ~ 15" for SCUBA at 850jum), current esti- 
mate softheJ}^£icalsiz£ofSM ra- 
dio ( Chapman et alJ 2004 Biggs & IvisonI 2008 ) and millime- 
ter (e.g. iTacconi et al T l2006r interferometric observations of 
small source samples. In particular. Biggs & Ivison (BI08 here- 
after) combined 1.4 GHz data from the VLA and MERLIN to 
produce high resolution radio maps of 12 SMGs detected in 
the Lockman Hole. The deconvolved sizes derived from these 
radio maps (obtained by fitti ng each map with an elliptical 
Gaussian) are consi stent with Chapman et alj (|2004|) (see also 
Muxlow et al ] l2005h and lTacconi et al. I (l2006l) . 



Figure [T] shows the effective radius and axis ratio for each 
of the 12 radio sources reported in BI08 (black circles). Note 
that, although the b/a interval [0.5, 1] used in optical lensing 
simulations contains 11 out of the 12 radio sources, there are 
several error bars that extend below its lower limit. In addition, 
the optical effective radius of 0.5" is not suitable to describe 
them. Since the source sample is too small (and the error bars 
rather large) to derive a reliable size distribution , we decided 
to use an effective radius close to the median of the sample. 
As a result, the size of the 1 .4 GHz radio emission produced by 
the SMG population has been characterized in the following by 
b/a randomly varying within [0.3, 1] and Rg = 0.25". 

Since SMGs se em to follow the FIR-radio correlation (e.g. 
Kovacs et al .120061) . their emission at both 1 .4 GHz and 850 yum 



is expected to be associated with massive star formation. This 
means that, as a first approximation, the same morphological 
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Fig. 2. Arc width probability distributions derived from a set of sources at Zs = 1 that were lensed by (the three projections of) the 
most massive cluster of the MareNostriim simulation at z = 0.3. The two panels illustrate the difference between selecting arcs 
with a length-to-width ratio d > do = 7.5 (left) and d > do - 10 (right). The black histograms were produced using the typical 
source size of SMGs at 1.4 GHz and 850 //m assumed in this work [Rs = 0.25" and b/a randomly varying in the interval [0.3, 1]). 
The blue histograms correspond to optical sources as usually described in lensing simulations {R^ - 0.5" and b/a € [0.5, 1]). 
The red dashed region shows the 150 mas resolution limit of the e-MERLIN radio interferometer at 1.4 GHz. 



parameters can be used to characterize the sizes and shapes of 
SMGs at radio and submm wavelengths. This choice of param- 
eters for the submm emission is also consistent with th e recent 
SMA observations presented by lYounger et al.l (l2008b . which 



this artificial increase of sources is corrected for by assigning 
a weight < 1 to each source, which is reduced at each new 
iterat i on step (see Miralda-Escude|[l993 : Bartelmann & Weiss 



have partially resolved the 890 yum emission of a SMG (GN20) 
for the first time (see Figure [U. 

If we wish to characterize gravitational arcs via their 
length-to-width ratio and make comparisons between obser- 
vations and theoretical predictions in an unbiased way, it is 
crucial to resolve their width. To address if the resolution pro- 
vided by radio and submm instruments could be an issue for 
arc statistics studies, we investigated the width distribution of 
arcs produced by a population of sources that is being lensed 
by a galaxy cluster. In particular, we used the most massive 
lensing clus ter at z - 0.3 fro m the MareNostrum cosmological 



simulation (IGottlober & Yepes.2007.) . a large «-body and gas- 
dynamical run, whose lensing properties recently have been 
studied by Fedeli et al. (2009, in preparation). The mass dis- 
tribution of this cluster was projected along three orthogonal 
directions, for which we der ived deflection angle maps b y stan- 
dard ray-tracing techniques ( Bartelmann & Weisslll994 ). Then, 
a set of sources at Zs = 1 (source redshift at which the lens- 
ing efficiency peaks for lenses at z ~ 0.3) with R^ - 0.25" 
and axis ratios randomly varying in the interval [0.3, 1] was 
lensed through the three projections. As usual in this proce- 
dure, the sources are preferentially placed near the lensing 
caustics following an iterative procedure to enhance the prob- 
ability of the production of large arcs. The bias introduced by 



1994[ iBartelmann et al.lll998l forfurflier details). 

The black histogram shown in Figure |2] corresponds to the 
width distribution of arcs derived from this simulation. For 
comparison, we have also included the corresponding distribu- 
tion derived from the parameters used in optical lensing sim- 
ulations (blue histogram). The two panels illustrate the dif- 
ference between selecting arcs characterized hy d > do with 
do - 7.5 (left) and do - 10 (right). As expected, reducing 
the source equivalent size produces a decrease in the width of 
lensed images. Note also that the behavior of the distributions 
is practically independent of the minimum length-to-width ra- 
tio used to select the arcs. The most important feature, however, 
is that both distributions drop to zero for widths below ~ 0.2", 
meaning that virtually no radio/submm (or optical) arcs have 
widths smaller than that value. At 1.4 GHz, ~ 0.2" resolu- 
tion is akeady accessible with MERLIN/e-MERLIN (~ 0.15"). 
Therefore, resolving the width of long and thin images for arc 
statistics studies is in principle already possible at radio wave- 
lengths. However, until the advent of ALMA, the ~ 0.75" reso- 
lution of the SMA at 950 /im will only be able to resolve a very 
small fraction of arcs produced by the most extended SMGs. 

Finally, we would like to stress two points related to the 
choice of source morphological parameters presented in this 
section. First, the most lumino us SMGs seem to be the result of 
merger processes (e.g.. lGreve et al.i,2005L .Tacconi et al 
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Fig. 3. Left panel. Redshift distribution of SMGs derived fro m the spectroscopic sample of IChapman et al.l (l2005h (cyan his- 
togram) and the photometric sample of lAretxaga et al.l (120071) (dark-grey histogram). The curves SPZ and PHZ correspond to 
the best fits provided by Eq. to the spectroscopic and photometric data, respectively. The curve CHM corresponds to the best 
Gaussian fit to the simple evolutionary model for SMGs used in CH05 (as quoted in CH05), normalized to the redshift interval 
[0, +oa\. Note that the binning adopted here is different from the one used in CH05, so the "redshift desert" in the redshift interval 
Z= 1 .2 - 1 .8 is not evident. However, changing the binning of the histogram changes little of the subsequent fit. Right panel. The 
cumulative distributions corresponding to the fits reported on the left panel, calculated by using Eq. Q. 



hence it is unlikely that their true shape is elliptical, as we as- 
sumed. However, if a merging source is lensed as an arc at a 
particular wavelength, irregularities in its shape and internal 
structure will not significantly change the global morphological 
properties of the arc, like the length-to-width ratio. What can 
happen is that the length-to-width ratio of an arc changes with 
wavelength because the emitting region of the source at those 
wavelengths have different sizes. For some of the wavelengths 
the source might even look like a group of small isolated emit- 
ting regions instead of a continuous one, which means that in 
the image plane it will be observed as a group of disconnected 
multiple images instead of a full arc. A very illustrative exam- 
ple of this scenario is SMM J04542-0301, an elongated region 
of submm emission which seems to be associated with a merger 
at z = 2.9 that is being lensed by the cluster MS045 1.6-0305 
(iBorys et al.l2003l:lBerciano Alba et al.ll2007l: Berciano Alba et 
al. 2009, in press). Until more complete information about the 
average structure of SMGs becomes available, we believe our 
approach to be the best that can be done. 

Second, since the radio sources studied by BIOS are brighter 
than 50 fjiy, the typical size derived from this sample might be 
different from the one that could be derived from fainter SMGs. 
Note, however, that the change of the cluster cross section with 
source size has a ve ry small slope for/?e between 0.2" and 1.5" 
dPedeli et al. I l2006h ." Therefore, deviations from - 0.25" 



(120061) . see Figure 6 of BIOS) will have a negligible effect on 
the derived number of arcs. 



3.2. Source redshift distribution 

A key point in trying to estimate the abundance of strong lens- 
ing features that are produced by the galaxy cluster population 
is the redshift distribution of background sources. Distributions 
peaked at higher redshift, or with a substantial high-z tail, will 
have in general more potential lenses at their disposal, and 
hence will produce larger arc abundances as compared to low- 
z-dominated distributions. In addition, the lensing efficiency for 
individual deflectors is also an increasing function of the source 
redshift. 

The most robust estimate of the redshift distribution of 
SMGs to date is based on the ~ 15 arcmin^ SCUBA survey 



within this interval (which is two times larger t han the inter- 
val that contains the sizes measured by BIOS and lTacconi et al 



carried out by IChapman et al.l (l2005h (CH05 hereafter). Radio 
observations were used to pinpoint the precise location of the 
submm detections, allowing the identification of optical coun- 
terparts that could provide precise spectroscopic redshifts. The 
final sample is composed of 73 SMGs with S50 i-im flux den- 
sities > 3 mJy and radio counterparts with flux densities at 
14 GHz>30/iJy. 

On the other han d, the SCUBA Half-Degree Extragalactic 
Survey (SHADES, iMortier et akl 120051: Ivan Kampen et al 
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2005 ) is the largest (720 arcmin^) 850 /im survey to datj*"! 
From their catalog of 120 SMGs, 69 have robust radio counter- 
parts with 5850 iim ^ 3 mJy and S 1.4 ghz ^ 20 ^ Jy. Photometric 
redshifts for this sub-sample were calculated bv lAretxaga et al 



( 12007b (AR07 hereafter) by fitting Spectral Energy Distribution 
(SED) templates to the available photometry at 850 fim and 
1 .4 GHz, including upper limits at 450 /urn (additional photom- 
etry at millimeter wavelengths was also used for 1 3 out of the 
69 sources). The histogram of the resulting photometric red- 
shift distribution, together with the spectroscopic one reported 
by CH05, are shown in Figure [3] The accuracy on the photo- 
metric redshifts derived by AR07 is Az ~ 0.65. Note that the re- 
quirement for a radio counterpart biases these two redshift dis- 
tributions against SMGs with z > 3, due to the less favorable K- 
correction in the radio compared with submm. Using a simple 
evolutionary model, CH05 estimated that the fraction of SMGs 
(S S5Q fim > 5 mJy) that is missing between z ~ 2.5 and z ~ 5 
in their spectroscopic redshift distribution is ~ 35%, a num- 
ber that is in agree ment with the fract i on of radio-unid e ntified 
SMGs reported in llvison et al] (l2002h : IChapman etaP (l2003h 
and the SHADES survey (AR07). The model evolves the local 
FIR luminosity function in l uminosity with inc reasing redshift 
following the prescription of lBlain et al.l (120021) . To account for 
the dust properties of SMGs, it also includes a range of SED 
templates th at have been tuned to fit their observed su bmm flux 
distribution (IChapman et al ] l2003l:lLewis et alJl2005h . 

In order to investigate the effect of this high-z tail on the 
predicted number of arcs, we used three different analytic ex- 
pressions to characterize the redshift distribution of SMGs in 
our calculations (see Figure [3]). One of them (CHM) corre- 
sponds to the best Gaussian fit to the distribution predicted by 
the CH05 evolutionary model, as quoted in CH05. The other 
two (SPZ and PHZ) were obtained by fitting the CH05 and 
AR07 histograms with the following anal ytic function, usua lly 
adopted for optical strong lensing studies ( Small et"al]|l995 ). 



P(Zs) = 



/3 



z; exp 







- 


.a) 







(4) 



where r{x) is the complete Euler-gamma function evaluated at 
X, Zq is a free parameter that broadly selects the position of the 
peak, and /? is another free parameter that defines the extension 
of the high-redshift tail. Note that, unlike the histograms, this 
function drops to zero at low redshift. However, the contribu- 
tion to the global lensing optical depth coming from sources 
at z < 1 is likely negligible, due to the geometric suppres- 
sion of lensing efficiency (see also the subsequent discussion in 
Section |5]). The resultant best-fitting parameters of these three 
functions are summarized in Table [1] 

The SPZ curve constitutes a good representation of CH05 
data not corrected for spectroscopic incompletenes0 The 
PHZ curve, on the other hand, does not describe the AR07 dis- 



The largest SMG survey to date (0.7 deg^) has being carried out 
at 1 . 1 mm with the AzTE C continuum camera mounted on the JCMT 
dAustermann et al.ll2009l) 

' ' Due to the lack of strong spectral features falling into the obser- 
vational windows, it is not possible to measure the spectroscopic z of 
sources in the range z = 1.2 - 1.8 (the so called "redshift desert") 



Table 1. Parameters of the redshift distributions presented in 
Figures |3] and m 



Nickname 


7"* 


^(1) 


rms^ ' 


<.p 


SPZ 


1.99 


2.81 




1.76 


PHZ 


2.51 


13.3 




2.18 


CHM 






1.30 


2.30 


SPZC'-" 
CHMC(*' 






0.64 
0.81 


2.05 
2.25 



best fit parameters of Eq. 

width of the best Gaussian fit. 

peak position of the distribution. 

the parameters associated with these 
distributions, provided by M. Swinbank, 
were misquoted in SW08 and CH05. 



tribution very well, failing to reproduce its z ~ 3 - 4.5 tail. This 
is a consequence of how the function in Eq. (|4]i is constructed. 
In particular, to accommodate the low-z part of the photometric 
histogram in Figure [3] the function needs to raise very steeply 
and hence, by construction, it must also drop steeply at high- 
z. Despite the function in Eq. ^ not being a good choice for 
fitting the photometric data, we nevertheless included the PHZ 
curve in our calculations because it highlights the consequence 
of choosing a distribution that is truncated at z ~ 3. Finally, 
the CHM curve allows us to predict the number of expected 
arcs if the CH05 histogram is corrected for spectroscopic in- 
completeness and high-z SMGs without radio counterparts (ra- 
dio incompleteness). We stress that at this stage we are not in- 
terested in using the best possible representation for the true 
source redshift distribution, but only to adopt a few motivated 
choices that broadly cover the range of realistic possibilities, 
in order to check the corresponding effect on the abundance of 
arcs. 

To show more clearly the difl'erent behavior in the high- 
redshift tail of our three choices, we present their cumulative 
distributions in the right panel of Figure |3] namely 



Piz) 



(5) 



In particular, when P(z) - 1 for SPZ, at z ~ 3.5, we still have 
P(z) - 0.8 for CHM, implying that ~ 20% of the SMGs still 
can be found at z > 3.5 using the latter distribution. 

In order to further show that this family of three functions 
cover all the reasonable possibilities, we have compared them 
with the predictions of one of the semi-analytical models that 
have been developed to explain the properties of SMGs (see 



Swinbank et al 



20081 and references therein). The histogram 
in the left panel of Figure |4] shows CH05 data coiTected for 
spectroscopic incompleteness. Note that the SPZ curve is very 
consistent with the semi-analytic model prediction (SWR), al- 
though both curves peak at slightly lower redshift (Azp - 0.3) 
than the best Gau ssian fit of the h i stogram (SPZC). However, 



as pointed out in ISwinbank et al.1 (120081) (SW08 hereafter) 



the CH05 distribution is expected to be uncertain by at least 
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Fig. 4. Left panel. Histogram of the redshift distribution of SMGs derived bv lChapman et al.l (120051) . corrected for spectroscopic 
incompleteness. This correction was implemented by interpolating the CH05 distribution in the region of the redshift desert 
(M. Swinbank, private communication). The SPZC line indicates the best Gaussian fit to the histogram. The blue solid line 
corresponds to the redshift distribution o f SMGs with 5850 /jm > 5 mJy and 5i.4ghz > 30 //Jy predicted by the semi-analytic 
model presented in lSwinbank et al.l (120081) . The SPZ curve presented in Figure[3]has being included for comparison. Right panel. 
The same histogram, also corrected for radio incompleteness using C H05 evolutionary mode l. The blue solid line corresponds to 
the redshift distribution of SMGs with 5 §50 nm > 5 mJy predicted in lSwinbank et al.l ( 120081) . The CHMC line indicates the best 
Gaussian fit to the histogram. The SPZ and CHM curves presented in Figure[3]have being included for comparison. 



Az ~ 0.25, which is the field-to-field variation between the 
seven sub-fields in the CH05 sample due to cosmic variance. 
Therefore, we can consider SPZ as a good representation of 
the current observations, despite the fact that it comes from a 
histogram that was not corrected for spectroscopic incomplete- 
ness. 

After the computations of the number of arcs were com- 
pleted, we became aware of the fact that the parameters quoted 
in CH05 for the best Gaussian fit to their simple evolutionary 
model (CHM, see Table [1]) were incorrect (M. Swinbank, pri- 
vate communication). As it is shown in the right panel of Figure 
m the CHM distribution has a higher-z tail as compared to the 
correct Gaussian fit (CHMC) and the prediction of the semi- 
analytical model (SWS). Since the true high-z tail of the red- 
shift distribution of SMGs is expected to be in between the 
cases considered in our calculations (SPZ, PHZ and CHM), 
and (as it will be discussed in Section 15.21 ) the final effect of 
the source redshift distribution on the number of arcs is small 
given the many uncertainties involved, we considered it unnec- 
essary to repeat the calculations for CHMC. 

3.3. Source number counts 

The final ingredient needed to estimate the number of arcs pro- 
duced by SMGs is the observed surface density of this source 
population, both at 1.4 GHz and 850 yum. Let no{S) be the dif- 



ferential number counts, defined as the surface density of un- 
lensed galaxies per unit flux density S . Integrating «o(5' ) over 
all fluxes above a given threshold, we obtain the respective cu- 
mulative number counts 



No(S) 



no{^)d^. 



(6) 



Let yu be the lensing-induced magnification of images on the 
lens plane, and ju+ 5 If P{fi+\do) is the magnification prob- 
ability distribution for sources that are imaged into arcs with 
d > do, then the magnified diff'erential numb er counts can be 
calculated as ( Bartelmann & Schneideij2001 ) 



«(5)= no - 
Jo 1/ 



S \ P{fi+\do) 



Hence, the magnified cumulative number counts read as 
Js Jo j"+ 



(7) 



(8) 



As can be seen in Eq. ([8]l, the lensing magnification bias has 
a twofold effect. On one side, sources that would be too faint 
to be detected without the action of lensing are amplified, and 
hence the respective images are brought above the detection 
threshold. On the other side, the unit solid angle is stretched 
by lensing magnification, implying that the number density of 
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sources is decreased. Which one of these two effects wins de- 
pends on the local slope of the unmagnified cumulative number 
counts. In particular, if N{){S) oc S^" and a > I, the number 
density of sources will be increased, while if a < 1 it will be 
decreased. It should be noted that, while for the unmagnified 
counts N^){S ) the flux density is derived by integrating the sur- 
face brightness over the area of the source, for the magnified 
counts N{S ) the integral is performed over the area of the re- 
sulting arc. 

Since our main motivation was to provide predictions for 
the abundance of giant arcs to be detected in surveys carried 
out with future instruments, we needed to provide the predicted 
number of arcs as a function of the surface brightness, instead 
of flux density. The reason is that we are working under the 
assumption that arcs are resolved structures, and therefore they 
are observed as extended objects. Under these circumstances, 
the flux integrated over the resolution element of the instru- 
ment (seeing, PSF, pixel or beam) is no longer the total flux 
of the source (as in the case of unresolved sources), and it 
may therefore be below the limiting flux although the arc as 
a whole is not. In other words, arc detectability under these cir- 
cumstances is not limited by the flux density but rather by the 
surface brightness. 

In order to take this into account, we had to convert the 
observed number counts as a function of flux density into num- 
ber counts as a function of surface brightness. Assuming that 
the size of sources is given by R^., and that the surface bright- 
ness is constant across it, then Nq{B) = No{S /ttRI). Note that, 
since the surface brightness is not affected by lensing, the mag- 
nification bias will manifest itself only through the solid angle 
stretching. Therefore, the cumulative magnified number counts 
(as a function of surface brightness) can be written as 



Table 2. Parameters for the 850 /urn differential number counts. 



N{B) = Nq{B) 



f 

Jo 



(9) 



Among other things, this implies that the magnification bias 
will always decrease the cumulative number counts, irrespec- 
tive of the shape of the unmagnified ones. 

In th e following we us ed the magnification distribution 
given by iFedeli et al.l (120081) . which is represented by the su- 
perposition of two Gaussians. In particular, we adopted the 
P(ju+\do) function for c/q - 10, but the result is virtually the 
same also for the case do - 7.5. Note however, that this (con- 
ditional) magnification distribution was computed for a back- 
ground population of sources that have different morphologies 
than SMGs (see Section 13.11 ). In principle, the bimodality of 
the magnification distribution is expected to b e preserved be - 
cause it only depends on the caustic structure ( iLi et alj|2005b . 
but it can be affected by the source morphology in two op- 
posite ways . On one hand, since SMGs are smaller than in 



Fedeli et al.l (120081) . we expect large arcs to form closer to the 
critical curves, and therefore to have larger magnifications on 
average. On the other hand, the fact that SMGs are more elon- 
gated will favor the formation of large arcs in regions of lower 
magnification. Given the uncertainties in other parts of the cal- 
culation, we consider that the use of a magnification distribu- 
tion derived for optical sources will have a marginal effect on 



Name 


S. (mJy) 


"o,* 


a 


/3 


DB('> 
SB*2) 
SM") 


4.30 ± 0.08 
4.22 


658 ± 48 
1039 ± 69 
970 


9 19+0.14 
^•^^-0.08 

-2.62 ±0.10 
-2.52 


6 22+°-5i 



best fit double power-law function, Eq. ilOb . 

best fit Schechter function, Eq. jilt . 

shallowest Schechter function consistent with the data. 

expressed in deg"- for DB and in deg"^ mJy"' for SB and SM. 



the derived number of arcs produced by SMGs. For a compre- 
hensive review of the many effects that could affect the estima- 
tion of arc abundan ces by galaxy clusters, see the discussion in 
(l2008h . 



FedeU et al 



3.3.1. Submm number counts 

For the latest and most complete estimate of the submm num- 
ber counts at 850 //m we refer to iKnudsen et al.l (l2008h (here- 
after KN08), who carried out a combined analysis of the counts 
derived from the Leiden SCUBA Lens Survey (LSLS) and the 
SHADES survey. With an area of 720 arcmin^ the SHADES 
survey is the largest blank-field submm survey completed to 
date, and therefore the least affected by cosmic variance. It 
provides the best constraints for t he submm number c ounts in 
the flux density range 2-15 mJy ( ICoppin et al.ll2006b . On the 
other hand, the LSLS survey targeted 12 galaxy cluster fields 
which cover a total area of 71.5 arcmin^ in the image plane. It 
provides the deepest constraints at the faint end of the submm 
counts (0.11 mJy, after correcting for the lensing magnifica- 
tion). In their analysis, KN08 used two functions to character- 
ize the combined differential number counts from both surveys: 
A double power-law. 



no{S ) = 



no,./5. 



and a Schechter function ( SchechteJll976l) . 



«o(5) = «o,J^| e-'"' 



(10) 



(11) 



Moreover, when fitting the observed cumulative number 
counts, they added the supplementary constraint that the inte- 
grated light well below 0.1 mJ y should not be higher than the 



extragalactic background light dPuget et al.llI996t iFixsen et al 
|l998V The resulting best-fit parameters are summarized in 
Tabled 

Figure |5] shows the (unmagnified) cumulative number 
counts derived from the best fit Schechter function (black solid 
line) and the best fit double power law function (red solid 
line) presented in KN08 using Eq. (|6]l. Note that, whereas 
both curves behave almost identically at the bright flux end, 
their predictions for the low flux number counts differ by a 
factor of ~ 2.5 at 0.1 mJy. Since the low flux tail of the 
submm counts dominates the number of SMGs that could po- 
tentially be lensed, we computed predictions for arcs produced 
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Fig. 6. The cumulative number counts of SMGs at 850 yum (top panels) and 1 .4 GHz (bottom panels), as a function of flux 
density (left panels) and surface brightness (right panels). The conversion from flux density into surface brightness was done 
assuming that the emission at both submm and radio wavelengths is homogeneous and have = 0.25". Thin lines correspond to 
S 1.4GHz = S 850 /im/50, while thick lines assume S i.4ghz = S §50 ;jm/100. The magnification pattern used to compute the magnified 
number counts is the one for arcs with d > do - 10. SM indicates the shallowest Schechter function that is consistent with the 
data, while DB represents the best fit double power-law (see the text for more details). 



by SMGs at 850 /im for the two following cases: (i) the shal- 
lowest Schechter function consistent with the combined LSLS 
and SHADES data (also shown in Figure |5]l and (ii) the best 
fit double power law function, hereafter refered to as SM and 
DB, respectively. The first one provides the minimum expected 
number of arcs consistent with observations, whereas the sec- 



ond one gives the number of arcs predicted by the best fit to the 
data (see Table|2]i. 

The cumulative number counts derived for these two cases 
as a function of flux density are shown in the top left panel 
of Figure |6] including the corresponding counts corrected for 
magnification bias using Eq. (O. In the same Figure, the top 
right panel shows the cumulative number counts as function 
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Fig. 5. Comparison between observed and predicted cumula- 
tive number counts. The red and black solid curves corresponds 
to the best-fit double powe r-law (DB) and t he bes t-fit Schechter 
function (SB) derived by iKnudsen et alj (l2008b for the com- 
bined 850 jum cumula tive number counts from the LSLS survey 



dKnudsen et al.ll2008h and the SHADES survey dCoppin et al . 
2006h .The black dashed line indicates the shallowest Schechter 



function consistent with the data of these two surveys (SM) 
The blue solid line indicates the cumulative nu mber counts pre- 



dicted by the semi-analytic model presented in lSwinbank et al 
(l2008h for SMGs with S sso /jm > 5 mJy. Model predictions for 



SMGs with 5 850 /jm > 5 mJy and S 1.4 ghz > 30 juJy, and SMGs 
with S&sQfjm > 5 mJy and 5'i.4ghz > 0.5 jiJy are indicated by 
the cyan solid and dashed lines, respectively. Note that the cyan 
dashed line and the blue solid line are almost indistinguishable. 



of surface brightness. The corresponding counts corrected for 
magnification bias (which will be used to compute the number 
of arcs) were derived using Eq. (|9]l. 

3.3.2. Radio number counts 

Figure|5]also shows the cumulative submm number counts pre- 
dicted by the SWOB model (blue solid line) compared with 
the results from diff'erent 850 yum SCUBA surveys. Note that, 
although the model tends to over-predict the counts at faint 
fluxes compared with the best fits provided by KN08 (red and 
black solid lines), it is consistent with the observational er- 
rors. The cyan solid line indicates the predicted counts for 
SMGs with radio counterparts assuming 5i4ghz > 30 yuJy. 
The fact that its shape is different from the shape of the blue 
solid curve is because current observations only detect radio 
emission from ~ 60% of the observed SMGs. However, if we 
allow the sensitivity threshold to go down to the //Jy level ex- 
pected for e-MERLlN, the SW08 model indicates that it would 



be possible to detect all the radio counterparts of SMGs with 
Sasojjm > 5 mJy (cyan dashed line). 

Since SMGs see m to follow the FlR/radio correlation (e.g. 
Kovacs et all 12006 ). the 1.4 GHz number counts of SMGs 



(which we need to predict the number of radio arcs) could be 
derived by scaling the 850 i-im number counts introduced in the 
previous section (DB and SM). As shown in Figure 7 of CH05, 
the ratio between the 850 lum flux density and the 1.4 GHz 
flux density shows a broad scatter (up to one order of magni- 
tude), which is probably a consequence of the strong influence 
of the dust temperature on the SEDs of SMGs. However, most 
of the points in this Figure are located between redshift 2 and 
3, and have an average 5850/jm/5i.4 ghz ratio between 50 and 
100. Therefore, we decided to use these two scaling factors to 
derive first order upper and lower limits of the radio number 
counts of SMGs. The resultant cumulative radio number counts 
are shown in the lower panels of Figure |6] Note that the values 
of 50 and 100 chosen for the submm/radio flux density ratio are 
meant to be indicative, since there are still many sources that 
display a ratio below 50 or above 100. The reader interested 
in results given by different values of this ratio can scale the 
curves appropriately in the upper panels of Figure |6] Also, ex- 
act numerical values can be made available by the authors upon 
request. 



4. Strong lensing optical depth 

To compute the total optical depth for lensed SMGs, we con- 
structed a synthetic cluster population composed of q = 1000 
cluster-sized dark-matter halos with masses uniformly dis- 
tributed in the interval [lO''', 2.5 x 10'^] Mq/i"' at z = 0. Note 
that it is not necessary to extract these masses according to 
the cluster mass function, since this is already taken into ac- 
count in Eq. ^ by weighting the cross sections with the func- 
tion n{M, z). The stru cture of each c luster is modeled using the 
NEW density profile dNavarro et al. 1995, 1996, 199?.) , which 
constitutes a good representation of average dark-matter halos 
over a wide range of masses, redshifts and cosmologies in nu- 
merical n-body simulations (,Dolag et al.ii2004) . Several studies 
of strong lensing and X-ray luminous clusters also show that 
these are well fitted by an NEW profile dSchmidt & Allenl2007l 



lOguri et al. 2009i) . This profile also has the adva ntage that its 



lensin g properties can be described analytically dBartelmann 



1996) 



To account for the asymmetries of real galaxy clusters, the 
halos are assumed to have elliptically distorted lensing poten- 
tials. However, instead of considering a single ellipticity value 
to describe all the synthetic cluster lenses, we derived an el- 
lipticity distribution from a set of n umerical clusters extracte d 
from the MareNostrum simulation dGottlober & Yepesll2007h . 
The strong lensing analysis required to generate this ellipticity 
distribution was taken from Fedeli et al. (2009, in preparation), 
as described in Section lTTI For each simulated cluster, the lens- 
ing analysis along three orthogonal projections was performed, 
computing the cross sections for arcs with d > do - 7.5 and 
sources with Zs = 2. For each of these projections, we found 
the ellipticity e of the NEW lens whose cross section is closest 
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Fig. 7. The distribution of NFW lens ellipticities fitting the 
cross sections of a sample of numerical clusters. The red dashed 
Une represents the best fit log-normal distribution, whose me- 
dian and dispersion are labeled in the top-right corner of the 
plot. 

to the cross section of the numerical cluster, i.e., we found the 
elUpticity that minimizes the quantity 



r(e) 



■ 0-7.5 



(12) 



where cr^,"^ is the cross section of the numerical lens, and 0-7 5(e) 
is that of the NFW lens for a given potential ellipticity e. 

Figure|2]shows the distribution of the ellipticities that min- 
imize the quantity r{e) in Eq. (fTSl i for all clusters in our nu- 
merical analysis. The dashed red line is derived by fitting the 
distribution with a log-normal function of the kind 



p(e)de — 



■ exp 



[ln(e) - In(eo)]' 



2cri 



dln(e), (13) 



where the best-fit parameters are eo = 0.31 and (Xe - 0.23. The 
ellipticity values used to characterize the potential of the syn- 
thetic NFW cluster lenses were then randomly extracted from 
the above distribution. 

Elliptical NFW profiles are a good representation of realis- 
tic cluster lenses only when the clusters do not undergo major 
merger events (jMeneghetti et al, 2003,). Since the merger ac- 
tivity of galaxy clusters is kno wn to have a sign i ficant effect 



on th e statistics of giant arcs (iTorri et alJ l2004t iFedeli et al 



20061) . it has to be taken into account in the construction of 
the synthetic cluster population. For thi s reason, we u s ed the 
excursion set formalisno developed by iLacev & Cold ( 1 1993b 



Also referred to as the "extended IPress & Schechtej i ll 9741) for- 
malism" 



(see also lBond et al.ll 19911 ISomerville & Kolattll 1999b to con- 
struct a backward merger tree for each model cluster at z = 
0, assuming that each me rger is binary (see discussion in 
Fedeli & Bartelmannl200"7l) . When a merger happens, the event 
is modeled assuming that the two merging halos (also described 
as elliptical NFW density profiles) approach each other at a 
constant speed. The duration of the me rger is set by the dy 



namical ti mescales of the two halos (see iFedeli & Barteknann 



(l2007n and IFedeli et al.l (l2008h for a detailed description of the 
modeling procedure). 

With the synthetic cluster population constructed in this 
way, the total average optical depth was derived by comput- 
ing individua l cross sections with the semi-analytic algorithm 
developed by IFedeli et al. especially designed to esti- 

mate the strong lensing cross sections of individual lenses in 
a fast and reliable way. The optical depth for a discrete set of 
lenses can be recast as 



2 nz^ ^ ^ nMi+[ 



dz, (14) 



where the masses M, ( 1 <ii^ q) have to be sorted from the low- 
est to the highest at each redshift step, and the quantity CTdoAz) 
is defined as 



o-doAz) = ^ [o-do{Mi,z) ■ 



■o-d,{Mi+uz)]. 



(15) 



This effectively means that, for all the clusters with mass be- 
tween Mi and M,+i, we assume the average cross section of the 
model d ark-matter halos w ith masses M, and M,+i. The algo- 
rithm of Fedeli et al. I (l2006l) for computing strong lensing cross 
sections consists of first assuming sources as point-like circles, 
and then introd ucing the effect of source ellipticity according to 
KeetonI (1200 lb . The source finite size is taken into account by 
convolving the local lensing properties over the typical source 
size. 

The total average optical depth is calculated by integrating 
Eq. ([T4l i over the source redshift distribution. Effectively, the 
pizs) weighting is avoided since we assigned individual source 
redshifts (randomly extracted from the distribution p{Zs)) at 
each of the q halos in the cluster sample and evolved their 
merger trees back in time until the respective source redshift. 
Given the large number of synthetic dark-matter halos used, 
this approach allows one to omit p{Zs) in Eq. dU when the red- 
shift integral is discretized. 

Despite the fact that the ellipticity distribution used to 
model the synthetic cluster population was derived for d > do - 
7.5, the best fit parameters of Eq. ([T3T l can be used to compute 
the cross sections for arcs also with c/o = 10 without compro- 
mising the results. The reason is that, although there might be 
a mild dependence of the distribution of lensing ellipticities on 
do, it is the overall caustic structure that defines the abundance 
of arcs above a certain do, regardless of its precise value. In 
fact, the criterion used to determine the ellipticity distribution is 
based on the similarity between the cross sections of the NFW 
lens and the numerical lens, which is an indirect way of com- 
paring the overall caustic structures produced by both kinds of 
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lenses. To verify this argument, the ellipticity distribution was 
re-computed using a criterion that is directly related to the caus- 
tic structure, that is, by defining the best fitting ellipticity as the 
one that minimizes the modified Hausdorff distanc^H between 
the critical lines of the numerical and NFW lenses. The ellip- 
ticity distribution obtained in this way is very similar to the one 
depicted in Figure|2](eo = 0.30, rms - 0.34). 

The median of the ellipticity distribution derived in this 
work using the lensing cross section (eo = 0.31, see Figure 



Table 3. Parameters of the arc redshift distributions shown in 
Figure ID 



is fu lly consistent with the one obtained by lMeneghetti et al 



(l2003h comparing the deflection angle maps (e ~ 0.3, using 
lenses placed at z ~ 0.3 and a source population at Zs = !)■ 
Even though these two criteria are arguably tightly related, the 
former is a quantity that is more directly related to observables 
than the latter, hence it is reassuring that they give compara- 
ble results. In addition, this result extends the previous one by 
quantifying the scatter around the median ellipticity and deal- 
ing with lenses that are distributed in redshift. 

5. Results 

5.1. Arc redshift distribution 

The redshift distribution of arcs with d > do (arc redshift dis- 
tribution hereafter) is expected to provide information about 
the redshift distribution of the background source population 
that is being lensed. However, it will also be distorted by the 
fact that the abundance of massive and compact galaxy clusters 
evolves with redshift, and that the lensing efficiency depends 
on the relative distances of sources and lens with respect to the 
observer. 

To assess the potential of this approach to gather informa- 
tion about the intrinsic redshift distribution of SMGs, we de- 
rived the arc redshift distribution associated with each of the 
three source redshift distributions used as inputs in our calcu- 
lations (see Figure |3]l. This was done by computing, for each 
input distribution, the average optical depth fd^izs) for several 
different values of Zs- That is, we excluded in the computation 
of the average optical depth those model clusters (and their 
respective chain of progenitors) whose associated source red- 
shifts were > z,. The resultant cumulative arc redshift distribu- 
tions, obtained after normalizing the optical depth for each Zs 
to the total average optical depth fda{+°°) - fdg, are shown in 
the right panel of Figure [8] The lines correspond to the best fit 
for each distribution provided by the simple function 



T"rf„(Zs) , / Zs 

= 1 - exp 



(16) 



where Zs,» indicates where the transition between the extrema 
and 1 occur, and y indicates how sharp this transition is. The 
best-fit parameters for each of the three input source redshift 
distributions are summarized in Table [3] The corresponding 
differential arc redshift distributions are shown in the left panel 
of Figure [8] 



This parameter constitutes one of the best ways to quan- 
tify the morphological differenc e between two set s of points 
jPubuisson & JaiiJ|l994) . See also iRzepecki et alj|2007l for a differ- 
ent application to gravitational lensing. 
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Once more, these arc redshift distributions have been 
shown for arcs with length-to- width ratio larger than do = 10 
only. As we verified, since the relative contribution of individ- 
ual model clusters to the total average optical depth is the same 
for both d() - 7.5 and do - 10, the resulting arc redshift dis- 
tribution also does not change significantly between the two 
choices. 

A comparison between Figure[8]and Figure[3]shows that, as 
expected, the arc redshift distributions reflect the general prop- 
erties of the source redshift distributions used as input, but there 
are also some noteworthy differences between them. For in- 
stance, the arc redshift distribution associated with CHM tends 
to zero at very low redshift, unlike in the case of the original 
CHM distribution. The reason is that low redshift sources do 
not produce many arcs, because (i) they have very few potential 
lenses at their disposal, and (ii) the lensing efficiency of those 
lenses is very low due to geometric suppression. This results in 
a lack of low redshift arcs in the distribution, which shifts its 
peak to higher redshifts compared with the CHM peak (from 
Zp = 2.3 to Zp ~ 3.2). Similarly, the peak of the arc redshift 
distribution corresponding to SPZ is shifted from Zp = 1 .76 to 
Zp >2. 

On the other hand, the arc redshift distribution correspond- 
ing to PHZ is not significantly shifted but visibly narrowed. As 
in the previous cases, low-redshift sources are removed from 
the distribution because they cannot produce arcs, but the dis- 
tribution could not shift at higher redshift because the input 
source redshift distribution is immediately truncated at Zs ^ 3. 
In other words, there is too little room between the drop due 
to lensing efficiency and the one due to the cutoff of the in- 
put distribution to allow a significant shift in its peak, and the 
only possible consequence for the distribution is to shrink and 
increase the peak height in order to preserve the normalization. 

In general, it is apparent that the differences between differ- 
ent source redshift distributions are somewhat enhanced when 
it comes to the arc redshift distribution. Therefore, assuming 
that redshift information is available for arcs, this approach can 
in principle be used to obtain some information about the gen- 
eral characteristics of the source redshift distribution, although 
it will probably not allow one to distinguish between redshift 
distributions that are very similar. 

5.2. Number of radio and submm arcs 

In this section we present and discuss the main results of this 
work: the predicted number of arcs produced by SMGs at radio 
and submm wavelengths. To that end, we computed the total 
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Fig. 8. Left panel. Differential redshift distributions of SMGs producing large gravitational arcs, corresponding to the cumulative 
distributions reported in the right panel. Right panel. The cumulative redshift distributions of SMGs producing radio arcs with 
d > do - \Q for input source redshift distributions SPZ , PHZ and CHM (see the text for more details). Lines are the best fit 
functions given in Eq. ( fTSI l and Table [3] according to the labels in the plot. 



average optical depth for each of the three source redshift distri- 
butions presented in Section[32](PHZ, SPZ and CHM), and for 
arcs with length-to-width ratio higher than both do - 7.5 and 
da - 10. These quantities were then multiplied by the magni- 
fied cumulative source number counts presented in Section 1331 
(SM and DB), to obtain the arc number counts as function of 
surface brightness. The results, extrapolated to the whole sky, 
are shown in Figure|9] A detailed list with the predicted number 
of submm and radio arcs for different sensitivities is presented 
in Tables |4] and |6] respectively. 

Note that the arc number counts given by the source redshift 
distributions SPZ and PHZ are almost indistinguishable on the 
scale of Figure |9] iiTespective of the length-to-width threshold 
c/o adopted. As expected, SPZ produces more large arcs than 
PHZ because it peaks at higher redshift, but only by a factor of 
~ 7%. The redshift distribution CHM, on the other hand, pro- 
duces more large arcs than the other two by a factor of ~ 2. 
Since, as mentioned in Section l372l the distributions CHM and 
PHZ can be considered as upper and lower limits to the true 
redshift distribution respectively, we can conclude that the un- 
certainty introduced by the redshift distribution in the predicted 
number of arcs is less than a factor of two. 

In terms of the length-to-width threshold, the number of 
arcs predicted for c/q = 7.5 is larger than for lio = 10 (as was 
also expected). However, the ratio between the number of arcs 
with d > do -1 .5 and c/ > c/q = 10 is not exactly equal to the 
ratio in the respective optical depths, since the magnification 
distributions for the two kinds of arcs are also different (see the 
discussion in iFedeli et al ..,20081) . 



Finally, when it comes to comparing the results from the 
two adopted source number counts (DB and SM), we see that 
the difference in the predicted number of arcs is negligible for 
submm surface brightness limits greater than 5 mJy arcsec"^. 
However, at Bgso/jm - 0.5 mJy arcsec"^, the function DB pre- 
dicts 2 times more arcs than SM, and the difference becomes a 
factor of 5 for 0. 1 mJy arcsec"^. In the radio domain, the differ- 
ence between DB and SM is negligible for Bn ghz — 50 ^Jy 
arcsec"^, a factor ~ 2 for 10 piy arcsec"^, a factor ~ 3 for 5 piy 
arcsec"^ and a factor ~ 8 for 1 piy arcsec"^. Therefore, the un- 
certainty in the predicted number of arcs is clearly dominated 
by the uncertainty of the source number counts at the faint sur- 
face brightness end. 

Considering an all-sky submm survey with enough resolu- 
tion to resolve individual arcs with d > do - 7.5 (~ 0.2"), 
and a sensitivity of Z?850 /jm - 1 mJy arcsec"^, our calculations 
predict between ~ 500 (PHZ SM) and ~ 600 (SPZ DB) arcs 
with a signal-to-noise ratio (SNR) larger than 5. In the case of 
d > do - 10, the expected number of arcs would be between 
200 and 250. If the sensitivity is reduced to 0.1 mJy arcsec"^^, 
these predictions can vary between ~ 3400 and ~ 8300 for 
do = 7.5, and between 1400 and 3600 for da = 10. 

In a similar way, an all sky radio survey with a sensitivity 
of Bi.4 GHz = 0.1 mJy arcsec"^ would detect between 8 and 50 
arcs for do - 1.5 (between none and 25 for do = 10) at SNR 
> 5, with the main uncertainty given by the Sg,5o /jm/S \aghz 
ratio used to obtain the radio number counts by scaling the ob- 
served submm counts. If the limiting radio surface brightness 
is reduced to ~ 10 piy arcsec"^, the predicted number of arcs 
could be increased to ~ 500 - 1300 for do = 7.5 (200 - 600 for 
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Fig. 9. The total number of arcs with d > do - 7.5 (left panels) and d > do - 10 (right panels) that are predicted to be observed in 
the whole sky above the surface brightness reported on the abscissa. Results for each of the three input source redshift distributions 
as well as both source number counts adopted in this work are shown, according to the labels. The two top panels refer to surface 
brightness at 850 fim, while the bottom panels refer to 1.4 GHz. In the bottom panels, thin lines refer to S 1.4 ghz = S^so //m/50, 
while thick lines refer to Siaghz = 5 §50 ^m/ 100. Please note the difference in scale on the abscissa between top and bottom 
panels. 



do — 10). In order to have excellent statistics with a few thou- 
sand giant arcs with SNR > 5, it would be necessary to go as 
deep as 5 i .4ghz = 1 j^/Jy arcsec"-. The largest number of arcs is 
given by the CHM source redshift distribution, which is about 
a factor of 2 larger than the number of arcs predicted by SPZ 
and PHZ. 



It is plausible that future radio and submillimeter surveys 
of galaxy clusters would focus on the most massive objects, 
since the center of attention of many multiwavelength studies 
is on X-ray bright clusters. To roughly evaluate the effect of 
this kind of selection, we re-computed the optical depths by 
including only those clusters in our synthetic population with 
mass M > 5 X IO'^^Mq/i-'. The resulting arc number counts 
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Fig. 10. The total number of arcs with d > do - 7.5 (left panels) and d > do - 10 (right panels) that are predicted to be observed 
in the whole sky above the surface brightness reported on the abscissa. Only clusters with mass are included 

in the calculations. Results for each of the three input source redshift distributions and both source number counts adopted in this 
work are shown, according to the labels. Top panels refer to submillimeter number counts, while bottom panels refer to radio 
number counts (notice the different scale on the horizontal axis). The difference between thin and thick lines is as in Figure|9] 



are presented in Figure [TO] where we used the same scale and 
line types as in Figure |9]to ease comparison. The correspond- 
ing numbers of submm and radio arcs predicted for different 
sensitivities are listed in Tables |5] and [T] respectively. The most 
interesting feature about these new plots is that the relative dif- 
ference in the arc number counts produced by the CHM dis- 
tribution on the one hand, and PHZ and SPZ on the other, is 
significantly reduced. This is due to the fact that massive clus- 



ters are found mainly at low redshift, hence lower-mass higher- 
redshift lenses, that are accessible only to the CHM distribu- 
tion, become unimportant. 

The order-of-magnitude reduction in the abundance of gi- 
ant arcs when focusing only on clusters with mass M > 5 x 
1O''*M0/!"' is in agreement with the fact that the bulk of the 
lensing signal actually comes from low-mass clusters, since the 
optical depth is obviously dominated by the lowest mass ob- 
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jects that are capable of producing a non-vanishing cross sec- 
tions (see Eq.[T]i. Note also that, unlike the previous case, PHZ 
produces slightly more arcs than SPZ. This is due to the fact 
that we are including in the calculations only low-z clusters, 
and PHZ actually has more sources with, e.g., Zs > 1 than SPZ 
(see Figure[3]). 

Considering again an all-sky submm survey with enough 
resolution to resolve individual arcs with do - 7.5 and sensi- 
tivity of Bgso/jm = 1 mJy arcsec"^, we predict ~ 20 arcs with 
SNR > 5 (8 if ^^o = 10). This number can increase to 100 - 250 
(40 - 90) if Bgso/im ^ 0.1 mJy arcsec"^. In the case of a radio 
survey, it would be necessary to push the limiting flux density 
down to ~ 1 yuJy arcsec"^ to detect few hundred arcs (between 
40-200ift/o = 10). 



6. Comparison with previous work 

The probability of strong gravitational lensing due to back- 
ground sources at radio and submm wavelengths has been a 
poorly studied issue in the past years. In addition, the few 
works available in the literature usually involve cosmological 
models, deflector mass ranges, and modeling approaches for 
the source and lens populations that are different from the ones 
used in the present work, making the comparison between them 
difficult and often not possible. With this note of caution, we 
tried however to make some of these tentative comparisons in 
the following. This required us to repeat the calculations in 
the last section using the number counts as a function of flux 
density instead of surface brightness. The results are shown in 
Figure [TT]only for arcs with d> do - 10. 



6.1. Submm wavelengths 

Given the poor resolution of current submm instruments, obser- 
vational studies of submm arcs have not been possible so far, 
with only one submm arc candidate reported until now. This arc 
is supposed to be the brightest region of the extended submm 
source SMM J04542 - 0301, located in the core of the clus- 
ter MS0451.6 - 0305 dBoryset all 120041: iBerciano Alba et al 

2009 



2007t Berciano Alba et al. 2009, in press). However, higher 
resolution submm observations are required in order to con- 
firm this hypothesis, and meaningful estimates of the length 
and width of this source, unfortunately, cannot yet be done. 
On the theoretical side, gravitational lensing of SMGs due to 
galaxy clusters was first studied by iBlain ( 1997]), using a cir- 
cularly symmetric model of the cluster A2218 and the submm 
counts predic ted by different galaxy evolution models. More re- 
cent studies (ICooray 1999; Paciga e t al...2009.) also have been 
focused on predicting the number of submm lensed sources, but 
predictions for the abundance of submm arcs have never been 
attempted before. 

(l2009l) . the au- 



For instance, in the work of iPaciga et al 



thors employ the strong le nsing analysis of the Millennium 



Simulation perfor med by iHilbert et al.l (l2007h (see also 



Hilbert et al.ll2008h in order to compute (i) the average mag- 
nification of SMGs as a function of flux density, and (ii) the 
contribution to the differential number counts given by sources 



with different redshifts and magnifications. Hence, their results 
cannot be compared with ours in a straightforward way. 

The only work w ith which we c ould try a tentative com- 
parison is the one by Cooray ( 19991) (C099 hereafter), where 
the author provides number counts of (among others) gravi- 
tationally lensed submm sources as a function of their mag- 
nification. The clusters were modeled as Singular Isothermal 
Sphere (SIS henceforth) density profiles, which means that the 
image magnification equals its length-to- width ratio, as long as 
sources are circular and point-like. Since the sources that we 
are using are not circular, nor point-like, the following com- 
parison should be taken with caution. The background submm 
sources were described by means of the redshift and num- 
ber distributions observed in the Hubble Deep Field (HDF). 
Using a ACDM cosmology, and considering cluster lenses with 
M > 8.8 X IO^^Mq/;-', C099 predicted ~ 500 submm sources 
in the whole sky with lensing magnification larger than 10 and 



S 



850 iim 



2 mJy. 



On the other hand, our results indicate that, if no mass se- 
lection is applied and for S'ggo^/m > 2 mJy, we should find 
a few thousand arcs with d > do - 10 in the whole sky. 
Restricting the cluster mass range to M > 5 x IO'^Mq/j"' re- 
duces the number of arcs to ~ 100 at most. Matching the mass 
range of C099 would only reduce the number of predicted arcs 
even further, hence being discrepant with C099 predictions. 
Assuming that the source population we are considering is the 
same, we ascribe at least part of this disagreement to the fact 
that C099 considers a very high normalization of the power 
spectrum (erg - 1.2), which was based on old stu dies of the 
cl uster temperature fu nction ( Viana & Liddlel[T996 ). As shown 



Fedeli et al.l (120081) . this results in greater strong lensing op- 



tical depth, and therefore to an over-prediction of the number 
of lensed images. 



6.2. Radio wavelengths 

The first obs ervational search for radi o arcs in galaxy clusters 
dates back to lBagchi & Kapahil ( ll995h . where the authors con- 
sidered a cluster sample of 46 objects with z < 0.3 observed, 
among others, at a wavelength of 20 cm with a sensitivity of 
~ 1 mJy. They claimed a systematic tangential alignment of 
radio images with respect to the cluster centers, and concluded 
that these are arclets or even giant arcs produced by background 
flat-spectrum radio sources. While their resolution was too poor 
to reliably measure arc morphologic al properties, the s e find - 



ings were subsequently questioned bv lAndernach et al 



who performed a similar analysis on the lAbell et al 



19981) 
1989h 



cluster sample (ACO henceforth) with sources taken from the 
FIRST (F aint Images of the R adio Sky at Twenty-centimeters) 
catalogue (iWhite et al ] ll997h . finding no evidence for a prefer- 
ential alignment of radio images in the core of massive galaxy 
clusters. They also report a rough estimate of the abundance 
of strong lensing events that should be seen in ACO clusters, 
finding that for a statistically significant detection, the limit- 
ing flux density should be lowered to < 0.1 mJy, in agreemen t 
with previous, more detailed estimates ( IWu & Hammeij[l993l) . 
More recentlv. Phillips et al..(.2001i) searched for strong lensing 
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Fig. 11. The number of arcs with d > do - 10 predicted to be observed in the whole sky as a function of the submm flux 
density at 850 //m (top two panels) and radio flux density at 1.4 GHz (bottom two). The left panels refer to all clusters in the 
sample, while the right ones shows results when only clusters with mass M > IO'^Mq/i ' are considered. Different line styles and 
colors refer to different number count models and source redshift distributions. In the two bottom panels, the thin lines assume 
S 1.4GHz = ^850 f,m/5Q, while the thick lines refer to S i.4GHz = ^850 ,jm/100. 



events with large angular separations in the FIRST catalog, but 
none of their candidates turned out to be a real gravitational 
lens. 

About 40% of the clusters included in the sample of 
Bagchi & Kapahil (1 1995b have velocity dispersion ~ 1000 



km s . Assuming that cluster galaxies have the same veloc- 



ity dispersion of th e dark-matt er particles (iGao et all l2004t 



In „ 

Biviano et al. 2006t Faltenbacher & Diemandll2006l see how- 



ever ICohn et al.llioool). and ad opting the simulation-calibrated 
scaling relation of lEvrard et al. (2008) . this corresponds to a 
mass of at least ~ 8 x IQ^'^M^h'^. The rest of their clusters 
should have a mass > 4 x IO'^^Mq/i"'. Assuming a limiting flux 
of ~ 1 mJy at 1.4 GHz, our predictions in Figure[TT]give < 100 
giant radio arcs in the whole sky when the entire cluster pop- 
ulation is considered, and only < 1 inside clusters with mass 
M > 5 X IO'^^Mq/i Therefore, we find it quite unlikely that 
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the detection claimed bv lBagchi & Kapahil (1 1995b is due to the 
radio emission from SMGs. 

The previous considerations suggest the need to go to flux 
densities fainter than 1 mJy in order to detect large radio arcs 
in galaxy cluster s. This kind of implication is also supported 
by the results of ICooravl ( 1 1998b . which showed that the radio 
emission corresponding to large optical arcs in three out of the 
four cases he studied is < 0.5 mJy. In the only optical arc with a 
secure radio counterpart in that work, namely arc AO in cluster 
A370, the radio counterpart does not seem to be arc-shaped, 
which may be due to the radio emitting region being offset with 
respect with the optical emission or to resolution issues. 

On the theoretical side, the statistics of radio arcs was 



first investigated by IWu & Hammed (Il993h (WU93 hereafter). 
T heir predictions were mad e using the evolutionary model 
of iDunlop & Peacockl (Il990l) to describe the radio luminosity 
function (dominated by starforming galaxies and AGNs prefer- 
entially located at z < 1) at 2.7 GHz assuming a cosmolog- 
ical model with Q.mfl = 1. Conversely, we focus on the ra- 
dio counterparts of SMGs (preferentially located at z ~ 2) at 
1 . 4 GHz in a more modern ACDM cosmology. As pointed out 



m 



Bartelmann et al.l (119981) . the optical depth for optical arcs 
produced by sources at Zs ~ 1 grows by about one order of 
magnitude in going from an Einstein-de Sitter universe to a 
ACDM one. Therefore, we can use this prescription for com- 
parison between both predictions. Note that, since the SED of 
SMGs is rathe r flat at radio wavelengths (spectral index ~ 0.7, 
Condon|[l992l) . we expect only a small change in the flux den- 
sity between 2.7 GHz and 1 .4 GHz. 

In their calculations, WU93 only considered clusters with 
z < 0.6 and a galaxy velocity dispersion > 800 km s which 
corresponds to a mass > 3.5 x lO'^^Mo/i"'. After applying the 
cosmology correction mentioned before, their predicted num- 
ber of giant radio arcs for the whole sky is (i) ~ 1000 if 
Sij GHz > 10 yuJy, (ii) < 200 if 5 2.7 ghz > 0.1 mJy and (iii) 
very small if 52.7 ghz > a few mJy (such that none should be 
detected in surveys in the literature at that time). Considering 
only clusters with mass > 5 x IO'^Mq/i"' (which have a red- 
shift range comparable to the one used in WU93), our predicted 
number of radio arcs with d > do - 10 produced by SMGs in 
the whole sky is (i) ~ 100 - 1000 if 5 1.4 ghz > 10 yuJy, (ii) few 
tens if 5 1.4 ghz > 0.1 mJy and (iii) very small if 5 1.4 ghz > a 
few mJy. 

Our results may thus seem compatible with those of WU99, 
considering that the number density of SMGs is certainly 
smaller than the number density of the entire radio source pop- 
ulation. This inference is however not conclusive, since WU99 
are using sources at rather low-redshift. We cannot say whether 
considering their same redshift distribution and number counts 
would lead to a discrepant result. 

In C099 there is also a study about strongly magnified ra- 
dio sources, analogous to the submm sources. It is found that 
~ 20 sources with amplification larger than 10 should be found 
with 5 1.4 GHz ^ 10 /iJy. For the same parameters, we find in 
our high-mass cluster study a number ranging from ~ 100 up 
to ~ 1000, which is much larger than the findings of C099. 
This discrepancy should be even more enhanced if the two cal- 
culations are reduced to the same a-%. The origin of this dis- 



crepancy is not clear, although it might be related to the higher 
mass threshold that they adopt, and the fact that the redshift 
distribution considered by C099 peaks at z ~ 1 instead of at 
z ~ 2. It should also be noted that C099 attribute their find- 
ing many fewer radio arcs than WU93 to the different number 
count evolution adopted. 

Our calculations of the abundance of radio and submm arcs 
are more accurate than the works discussed above for different 
reasons, mainly the different modeling of the cluster popula- 
tion. To start with, in the literature, investigators often consider 
all le nses as isolated, spheri cally symmetric density distribu- 



tions. |Wu_&jlannne2 (ll993|) investigate the effect of elliptical 
mass distributions, but only on the magnification pattern and 
not on the efficiency for the production of large arcs. On the 
other hand, we included the effect of asymmetries, substruc- 
tures and cluster mergers, that all have been found to be impor- 
tant to augment arc statistics. Next, we used an NEW density 
profile to model individual lenses, which is a good representa- 
tion of average dark-matter dominated objects like galaxy clus- 
ters, while other works have often considered SIS or SIS-like 
profiles, which are more suitable for galaxy lensing. While for 
a SIS lens model the image magnification equals the length-to- 
width ratio, it is known to produce fewer grav itational arcs with 



respe ct to the more realistic NEW profile (IMeneghetti et al 
2003h. 



7. Summary and conclusions 

The advent of the high resolution submm facilities ALMA and 
CCAT, and the major technological development that radio in- 
terferometry is currently undergoing (e.g., e-MERLIN, EVLA 
and SKA) will make possible the study of radio and submm 
giant arcs produced by clusters of galaxies. In particular, the 
study of giant arcs produced by submm galaxies (SMGs) seems 
particularly promising for at least two reasons. 

• It provides the opportunity to detect and spatially resolve 
the morphologies and internal dynamics of this population 
of dust obscured high-redshift star-forming galaxies, which 
is very difficult to study in the optical. 

• It can provide information about the formation and evolu- 
tion of the high redshift cluster population, by means of arc 
statistics studies. 

To assess the prospects for these kind of studies, we pro- 
vided theoretical predictions on the abundance of gravitational 
arcs produced by the SMG population at radio ( 1 .4 GHz) and 
submm (850 jim) wavelengths, greatly improving the accu- 
racy of the results with respect to the first few studies car- 
ried out a decade (or more) ago. The advantage of radio ob- 
servations is that the angular resolution and sensitivity pro- 
vided by interferometers like e-MERLIN and EVLA are al- 
ready (or will very soon be) at the level required for these kind 
of studies. However, these frequencies do not benefit from the 
same favorable K-correction as mm/submm wavelengths do, 
which will make the latter a more interesting tool for studying 
high-redshift clusters as soon as the resolution of sub-mm ob- 
servations is good enough. The calculation of the number of 
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arcs produced by a background source population requires four 
main ingredients: (i) the source shape and size (ii), the source 
redshift distribution, (iii) the cumulative source number counts, 
and (iv) a model of the cluster population. 

The model of the cluster population used in this work was 
based on the extended lPress & Schechten ( 119741) formalism and 
made use of an NFW dentity profile to describe each cluster 
lens. It also included the effect of asymmetries, substructures 
and cluster mergers, which have been found to play an impor- 
tant role in arc lensing statistics. 

Based on current radio/CO observations and the FIR/radio 
correlation, we have characterized the typical size of the radio 
and submm emitting regions of SMGs with an effective radius 

- 0.25", and an axis ratio that varies within the interval 
b/a e [0.3, 1]. Resolving all the arcs produced by this kind of 
sources will require ~ 0.2" resolution. 

Since the most accurate redshift distribution of SMGs avail- 



able dChapman et alj|2005h is based on observations of the ra- 
dio detected members (which is biased against z > 3 sources), 
we used three different functions to quantify the effect of a high 
redshift tail on the predicted number of arcs. The results indi- 
cate that this effect is less than a factor two if we consider all 
simulated clusters during the calculations, and negligible if we 
only consider massive clusters. 

The submm source number counts used in this work cor- 
respond to the joint fit of the (bright) SHADES survey counts 
and the (faint) Leiden SCUBA Lens Survey counts presented in 
Knudsen et al.l ( l2008h . To account for the uncertainty in the low 
flux end, predictions were made for the best fit to the data, and 
the shallowest fit consistent with the data. Note that, although 
only ~ 60% of the observed SMGs have being detected in ra- 
dio, the next generation of radio interferometers will be able to 
detect the radio counterparts of all SMGs with S 85o ^ 5 mJy. 
Therefore, the radio number counts of SMGs have been derived 
by scaling the submm counts using representative upper and 
lower limi ts of the S s'in am/S ] a r. H?, ratio for SMGs at z ~ 2 

(l2005h . 



taken from lChapman et al. 

Our calculations show that a submm all-sky survey with a 
sensitivity of 1 mJy arcsec"^ will detect hundreds of arcs with a 
5cr significance. In the radio, this number can be achieved with 
a sensitivity of 10 - 20 fiJy arcsec"^. Obtaining a statistically 
significant sample of thousands of arcs would require sensitivi- 
ties of 0.1 mJy arcsec"' in the submm and 1 juJy arcsec"' in the 
radio. However, if only massive clusters (M > 5 x IO'^Mq/i-') 
are considered in the calculations, the predicted number of arcs 
is reduced by about an order of magnitude. 

Besides the many uncertainties involved in the theoretical 
predictions presented here, the main challenge in designing a 
future survey for radio/submm arc statistics studies will be in 
finding the best compromise between survey area, depth and 
resolution, three issues that affect the arc detectability in differ- 
ent manners. We believe that this work provides a significant 
step forward in this direction. 
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Table 4. Predicted number of submm (850 /urn) arcs produced by SMGs for an all-sky survey, using all the synthetic cluster 
population to compute the optical depths (see Figure|9]l. 
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1550 
1442 
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4200 
3900 
475 



surface brightness limit used to determine the number of arcs, 
arc length-to-width ratio thresholds. 

redshift distribution and number count function, as introduced in Figures[3]and[6] 



Table 5. Predicted number of submm (850 jim) arcs produced by SMGs for an all-sky survey, using only clusters with M > 
5 X lO''* Mq/i"' to compute the optical depths (see Figure [TOli. 



850 iim 
^(2) 


5 mJy arcsec ' 
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17 8 
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75 25 
75 25 
83 33 


117 42 
125 42 
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308 108 
333 108 
383 133 



surface brightness limit used to determine the number of arcs, 
arc length-to-width ratio thresholds. 

redshift distribution and number count function, as introduced in Figures[3]and[6] 
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Table 6. Predicted number of radio ( 1 .4 GHz) arcs produced by SMGs for an all-sky survey, using all the synthetic cluster population to compute the optical depths (see Figure 
|9ll. Values with and without brackets correspond to S&so fjm/S 1.4 ghz = 100 and 50, respectively. 
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surface brightness limit used to determine the number of arcs, 
arc length-to-width ratio thresholds. 

redshift distribution and number count function, as introduced in Figures[3]and[6] 



Table 7. Predicted number of radio (1.4 GHz) arcs produced by SMGs for an all-sky survey, using only clusters with M > 5 x 10'"* MqIi ' to compute the optical depths (see 
Figure [Toll. Values with and without brackets correspond to Sgso fim/S 1.4 ghz = 100 and 50, respectively. 
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surface brightness limit used to determine the number of arcs, 
arc length-to- width ratio thresholds. 

redshift distribution and number count function, as introduced in Figures[3]and|6] 



